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A phenomenological approach is presented that allows one to model, and thereby inter- 
pret, photoemission spectra of strongly correlated electron systems. A simple analytical 
formula for the self-energy is proposed. This self-energy describes both coherent and 
incoherent parts of the spectrum (quasiparticle and Hubbard peaks, respectively). Free 
parameters in the expression are determined by fitting the density of states to experi- 
mental photoemission data. An explicit fitting is presented for the Lai-^Sr^TiOs system 
with 0.08 < x < 0.38. In general, our phenomenological approach provides information 
on the effective mass, the Hubbard interaction, and the spectral weight distribution in 
different parts of the spectrum. Limitations of this approach are also discussed. 



1. Introduction 

Photoemission experiments provide important information about the electronic 
single-particle excitation spectrum of solidsB For weakly correlated materials this is 
essentially given by the density of states (DOS) obtained by, e.g., deiisity functional 
theory in combination with the local density approximation (LDA).cl In many cases 
the agreement between LDA and experiment turns out to be very good. However, 
there is a class of materials wherjs the discrepancy between the measured and cal- 
culated spectra is significant .ulj'ffll3'Bflll30 For these strongly correlated electron 
systems there is a clear demand for new theoretical and computational approaches. 
The recently developed LDA+DMFT method, a combination of the LDA and the 



dynamica 
respect. IT" 



mean field theory (DMFT), has proved to be very successful in this 
100000 The LDA+DMFT method supplements the LDA by lo- 
cal correlations between d- or /-electrons .00 In the simplest case, namely, in the 
absence of long-range order and when the c orrelated bands at the Fermi level are 
sufficiently separated from other bandsjijlij the DOS of a correlated system is well 
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represented by the integral (Hilbert transform) 

Plda+dmftM = Im / cko ; , (1) 

7T J us — cj' — Ldmft(^) + «0 + 

where plda(w) is the DOS from the LDA calculation, and Sdmft(w) is the lo- 
cal self-energy calculated self-consistently within t^DJVIET^hepifi which includes 
correlation effects missing in the LDA approach.llail3allail3't3ll3 Non-local con- 
tributions to the self-energy cannot yet be implemented in this scheme. This will 
become possible in extensions of the DMFT, e.g., in th£L Dynamical Cluster Ap- 
proximation (DCA) and related computational schemesJia 

The aim of this paper is to describe another approach, phenomenological in 
nature, to model photoemission spectra of strongly correlated electrons. The mo- 
tivation is the following: the LDA+DMFT method is microscopic but requires an 
extensive numerical effort to calculate Plda(w) and Sdmft(^)- On the other hand, 
an analysis of various models of strongly correlated electrons within the DMFT 
has shown that certain features of the self-energy do not depend on the details 
of the model. Fermi liquid behaviour is seen in, e.g., Numerical Renormalization 
Group (NRG) calculations for the Hubbard model both at and away from half-filling 
(see Fig. (J2 ) below) and for the Periodic Anderson Model in the heavy-fermion 
regimeB3'E!J In these systems, the imaginary part of the self-energy consists of the 
w 2 -dependence for u> — > and T — > and two pronounced peaks at finite u>. On 
increasing the temperature, Fermi liquid behaviour can be rapidly destroyed, in 
particular close to the Mott transition,E3 or in general for systems with a very high 
effective mass. The imaginary part of the self-energy then goes over to a single 
and very broad peak centered approximately at the Fermi level (see, e.g., Fig. 5 
in Ref.[22]). This model- independence of the self-energy suggests the use of a uni- 
versal form for £(w) which depends only on a small number of phenomenological 
parameters. This replaces Edmft(w) in Eq. ([!]). 

Although the proposed scheme is phenomenological (the parameters in the self- 
energy being determined by fitting to the experimental data) we believe it to be 
useful for the qualitative interpretation and understanding of the experimental re- 
sults. The phenomenological self-energy obtained in this way can be used to 
deduce other quantities for the specific material, such as a linear specific heat coeffi- 
cient and the dynamical conductivity (under the assumption that vertex corrections 
are negligible). This approach would then serve as a unifying phenomenological de- 
scription of a variety of experimental results. 

Conceptually, such an approach is not new. It was used previously to fit and in- 
terpret, for example, the integrated photoemission data for Cai_ K Sr x VO_3,B and the 
angular resolved photoemission data in prototype Fermi liquid metalsEjca and high- 
temperature superconductors .El However, only the quasiparticle peak was fitted in 
these approaches. Here we provide an analytical expression for the self-energy which 
is appropriate for fitting the whole spectrum of correlated d- or /-electrons, where 
the Hubbard subbands and the quasiparticle resonance are present simultaneously. 
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Our approach is based on a sum of Lorentz functions. Matho has proposed another 
route based on the multi-pole expansion of the phenomenological self-energy.c3 It 
turns out that these two approaches are mathematically equivalent. In our approach 
it is possible to present a simple physical motivation for the form of the self-energy. 

Our phenomenological approach encounters certain difficulties which should be 
mentioned here. The most serious problem is connected with the description of 
multi-band systems. The strong electronic correlations originate from the localized 
nature of /- or d-orbitals. Hence, several bands might cross the Fermi level or 
be very close to it even when they are split by, e.g., a crystal field. In V2O3, for 
instance, the splitting between the ti g and e g bands is rather small. In such a case, 
each band which lies in the vicinity of the Fermi level would require a separate 
self-energy, which makes the number of fitting parameters twice or three times 
larger. Since the photoemission results are not orbitally resolved, an unambiguous 
fitting cannot be guaranteed in these cases. Without additional experimental input, 
a phenomenological approach for these cases is not adequate. 

In our paper we therefore concentrate on the experimental data for Lai-^Sr^TiOa , 
a system with degenerate ti g bands (see Sec. 4). Before that — in Sees. 2 and 3 — 
the phenomenological expression for E(cj) is introduced. The results of our paper 
are summarized in Sec. 5. 



2. Self-Energy 



We start with a heuristic derivation of the retarded self-energy for strongly corre- 
lated electrons (e.g., d-electrons) in the metallic phase which form a Fermi liquid 
state at low energies and temperatures. The DOS p(u>), calculated with this self- 
energy should consist of three parts: two wide incoherent parts (upper and lower 
Hubbard bands) and a coherent peak at or close to the Fermi level /i. 

The construction of a suitable self-energy expression is guided by the following 
idea. Let us start with a model for the spectral function j4 mo d(w) which is a sum 
of three Lorentz curves 



AnodM = 



Q 
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One peak (the quasiparticle peak) is centered at to = luq with spectral weight 
Q > and width 7 > 0. The other two peaks (upper and lower Hubbard peaks) 
are centered at w = ±1/2 and their widths are assumed to be both equal to T > 
(see Fig. (0)). The total weight of these two peaks is 1 — Q with the relative weights 
q/2 > and 1 — q/2 > respectively. A moc j(w) is normalized to unity. Apart 
from assuming the same width for both upper and lower Hubbard peaks, Eq. (||) 
is the most general sum of three Lorentz peaks. Since we have not yet specified 
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the position of the chemical potential fi, Eq. (^|) describes both symmetric and 
asymmetric cases. 




Fig. 1. Schematic plot of the three peaks of the model spectral function A mot j(<^). 

There is, of course, some arbitrariness in the choice of Eq. (||) for A mo a(uj). 
Lorentzians are used here to obtain simple analytical expressions. One can try to 
use, for example, Gauss or semi-elliptic forms of the DOS as well. However, in these 
cases either cannot be expressed analytically or is not a smooth function of 

U). 

The retarded local Green function <7 mo d(^) corresponding to the spectral func- 
tion (0) is 
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We also introduce a Green function 

SmodH = 7T-, (4) 

oj — 0Ji + zA 

which corresponds to a spectral function A^ od (o;) with a single peak centered at 
uj = lo\ and with the width A > 0. Note that all the above quantities (i.e. those 
with an index 'mod') do not correspond to any physical quantity but are introduced 
for the construction of a suitable self-energy only. 

Using the Green functions (||) and (Q) together with the Dyson equation for the 
self-energy T, mod (oj) = g" lod (w) _1 - 5 mo dM _1 we find 

-l 



J mod 



(iS) = UJ — 0J\ + iA — 



(1-0) 



w + »T-(l-g)| 



(5) 



This self-energy contains 8 parameters. The number of parameters may be reduced 
by imposing additional conditions which are discussed below. 

In order to preserve the Fermi liquid properties at low energy _we have to 
supplement the self-energy (||) by the condition ImS(u; = /z) = 0.E3 Then, as 
u> — > 0, ReE(oj) ~ —uj, and ImE(w) ~ — to 2 . In the high energy limit oj — > oo, 
ReS(w) ~ but ImX!(w) ~ a — with a constant a > 0, which means that 
the imaginary part of the self-energy changes sign. We have checked that 

this artefact becomes important only if / is much smaller than 7. However, in this 
weakly correlated limit plda usually reproduces the experimental data reliably and 
the corrections due to E(u>) are not necessary. In the strongly correlated limit the 
change of the sign in ImS(w) appears at high energies. We have therefore introduce 
a cut-off setting ImS(w) = whenever the self-energy becomes positive. 

The self-energy (||) is temperature independent. As noted in the introduction, 
the self-energy develops a peak at the chemical potential fj, for finite temperature, in 
particular close to the Mott transition. This effect is described phenomenologically 
by introducing a scattering part 

c 

SscattM = — — , (6) 

W — H + l^ls 

with two fitting parameters s and j s . The scattering part is not used in this paper 
but it might be important in systems far away from the Fermi liquid regime. 
Hence, the phenomenological self-energy takes the form 

Sfit(w) = S mod (w) + £ S cattM • (7) 

Before applying Eq. (]?]) to model experimental data, we check whether this 
form of the self-energy can reproduce the self-energies obtained numerically from 
the DMFT equations of the Hubbard model 

H =-t ( c h c 3<? + c ]a c i<?) + U Y1 c h c ^4i c H> ( 8 ) 

<ij>a i 

where t is the hopping matrix element between nearest neighbour sites and U is 
the local interaction energy for the electrons with antiparallel spins a. Figure H 



1, 2008 20:44 WSPC/Guidelines photo'070502 



6 K. Byczuk, R. Bulla, R. Claessen, D. Vollhardt 

shows the result of the fit to the spectral function, calculated by NRG with the 
microscopic parameters U = 4, n = — 1.4 and T = 0. We imposed the Fermi liquid 
condition Efit(w = /i) = 0. The scattering part S scat t(a;) was set to zero. The 
parameters determined from the fit-procedure are ojq = 0.02, 8 = 0.005, 7 = 0.005, 
T = 0.51, Q = 0.3, q = 0.72, and / = 4.5. We used the same bare DOS as in the 
NRG calculation, i.e. a semielliptic DOS with the width W = 2. This width sets 
the energy units in this fitting. The comparison shows both the possibilities and 
limitations of our phenomenological approach. Although the three-peak structure 
and the overall distribution of the spectral weight are described correctly, there are 
significant deviations regarding the form of the peaks. The main reason for this is 
that the peaks in the numerical data are not Lorentzian (see, e.g., the discussion of 
the form of the Kondo resonance in the single impurity Anderson model in [28]). 
The fit-procedure therefore cannot recover the dip at ui « 1.5, and compensates this 
by underestimating the width of the upper Hubbard peak. Also, the band filling 
determined from the fitted spectral function is larger by about 7% as compared to 
the NRG result. 

The same holds for the structures obtained for the self-energy as is seen in 
the lower panel of Fig. @. These are not, according to the NRG result, given by 
Lorentz peaks. Nevertheless, the general structure, i.e., the Fermi- liquid behavior 
for small frequencies and the two peaks at higher frequencies (oj — /i ~ —0.5 and 2.8) 
is reproduced correctly. The relative difference in the weight of these latter peaks 
is due to the particle-hole asymmetry in the parameters used for this particular 
calculation. 



3. Modeling Photoemission Spectra 

The phenomenological form for the self-energy Eq.(^|) discussed in the previous 
section can now be inserted in the Hilbert transformation for the DOS 

7T J LU - LU' - Sfit(w) + l0+ 

The direct (^direct) and inverse (inverse) photoemission intensities, within a con- 
stant transfer matrix approximation, are given by 

SdirectM - S J dw'R a {u - U')f (^jA) flit («>')> ( 10 ) 



-'inverse 



Pfit(w'), (11) 



where f[(co — [i)/k^T] is the Fermi-Dirac function with the chemical potential fi 
at the temperature k^T (in energy units) and R a (u>) = exp(— w 2 /2<t 2 )/v / 27tct is 
the apparatus function with the resolution a. Note that the resolution in direct 
photoemission experiments is typically about one or two orders of magnitude better 
than in inverse photoemission experiments. Sq and S' are constant prefactors. 
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Fig. 2. Upper panel: The fitted spectral function using the self-energy (^) (solid line) and the 
numerical spectra from NRG (dashed line). Lower panel: The phenomenological self-energy cor- 
responding to the fitting in the upper panel (solid line) and the self-energy from NRG (dashed 
line). Scales on the horizontal axis are the same. 

With good quality data for the direct and the inverse photoemission spectra on 
the same sample under the same conditions one can now determine the phenomeno- 
logical parameters in the self-energy (fjj) and determine, for example, the effective 
mass, the magnitude of the Hubbard interaction and the full frequency dependence 
of the self-energy. 

4. Example of Fitting 

We now turn to the phenomenological modeling of experimental photoemission 
data, as exemplified in the case of Lai^Sr^TiOs. In this compound, the 3d 1 elec- 
trons occupy degenerate t2 9 orbitals for which the crystal and Jahn- Teller splittings 
are very small. Moreover, the ti g band is well separated from the e g and p-bands. 
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Our phenomenological modeling can therefore be restricted to a single degenerate 
band. 

The available direct photoemission spectral! show the typical features of a strongly 
correlated metal, i.e., a quasiparticle peak and a lower Hubbard band. At low tem- 
perature Lai-^Sr^TiOs is an antiferromagnetic insulator for x = and a band 
insulator with an empty <i-band for x = 1. Consequently, the filling of the three- 
fold degenerate d-band decreases from n — 1 to n — in going from x = to x = 1. 
The antiferromagnetic insulator is stable up to x « 0.05 for T < Tm (?w = 112 K 
for x = 0), and for 0.05 < x < 0.08 an antiferromagnetic metallic phase appears 
with decreasing Tjv with increasing x. 

In the photoemission measurement, the quasiparticle peak is clearly visible in 
Fig. (|) for all paramagnetic samples with x = 0.08, 0.18, 0.28, 0.35.1 This co- 
herent peak is suppressed in the antiferromagnetic metallic phase (x = 0.06), and 
vanishes in the insulating regime (x = 0.04). The wide incoherent peak — the lower 
Hubbard band — is present for all x-values. These features, in particular the lower 
Hubbard band, cannot be explained within the LDA approacho 

Here we have used the experimental data! to determine the phenomenological 
parameters in the self-energy (Q). However, since we only know the occupied part 
of the spectrum it is impossible to determine unambiguously the absolute value 
of /, corresponding to the distance between the lower (occupied) and the upper 
(unoccupied) Hubbard bands. Only a relative value 1/2 + /i can be determined. In 
order to make the problem tractable we have to reduce the number of parameters. 

To this end we set / = 5 eV since this is the value found in other theoretical 
studies. ta Also we assume that the t2 g band is three-fold degenerate and that x = 1 
corresponds to 1 /3 filling. With decreasing x the filling of this band is lowered and 
is found to be n = (1 — x)/3, after normalizing the DOS to unity. This filling for 
a given x is used as another constraint on the parameters in our self-energy (Q). 
In other words, the parameters are adjusted such as to obtain the correct filling 
n = (1 — x)/3 of the d-band. We do not use the Fermi liquid constraint because the 
experiment was performed at finite temperature and close to the Mott transition, 
so that deviations from ImS(w = 0) = are expected to be significant. 

Within these assumptions the best fit is obtained by minimizing the mean square 
deviation between the theoretical and experimental data.ta The apparatus resolution 
function is taken to be a Gauss function with the variance a — 0.035 eV and the 
temperature is T = 20 K. The bare DOS is used as for LaTiOa from Ref. [29]. One 
example is shown in the lower panel of Fig. (|§|) . 

For comparison, we have also fitted the coherent part of the spectrum using the 
phenomenological self-energyd 

„ / n auj b 

SqpH=5— " —7- 12 

oj + la uo + ib 

This form of £(w) is useful to fit the quasiparticle peak but does not describe the 
Hubbard bands. In situations with no clear separation of these two structures, a 
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Fig. 3. Upper panel: Four curves from the photoemission experiments^) for Lai—^Sr^TiOa. Black 
dots show digitized points and small wiggles ^ase due to digitization errors. The inset shows the 
DOS from LDA calculations for this systemO Lower panel: Example of fitting the self-energy 
([?]) (solid line) and the self-energy ( |l^ ) (dot-dashed line) to the experimental result (dashed thick 
line) with x = 0.08. 



fit using Eq.(|12|) obviously involves some arbitrariness. This is in contrast to our 
fit formula for the self-energy which does not require the two structures to be well 
separated. 

As we can see in Fig. (||) both self-energies (0) and ( |l2l) provide a similar descrip- 
tion of the quasiparticle peak. However, the self-energy (Q) also fits the high-energy 
feature of the spectrum. In Fig. (Q) we show the total spectral functions (upper 
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-4.0 -2.0 0.0 2.0 4.0 6.0 8.0 

CO 

Fig. 4. Upper panel: The spectral functions for different x obtained from fitting to experimental 
data, as described in the text. Lower panel: Imaginary part of the self-energy for different x. The 
inset shows the behavior of l/Z, extracted from the real part of the self-energy, as a function of 
x. 

panel) for four different values of x obtained with this procedure. The lower panel 
presents the imaginary parts of the self-energy (||). As we can see, with decreasing 
x the number of states below the Fermi level is reduced and more spectral weight 
is pushed into the upper Hubbard peak. Furthermore, the quasiparticle peak is 
suppressed for higher x. These trends are also mirrored in the behavior of the self- 
energy for different x. The absolute value of ImE at the Fermi level increases with 
x. This increase might be attributed to the enhanced scattering due to the ran- 
domness introduced by the Sr atoms in the system. We also calculated the Z-factor 
{1/Z = 1 — 9RcS(cj = ii)/du>) which in Fermi liquid theory is related to the ef- 
fective mass of the quasiparticles {1/Z ~ m* /m). In the inset to Fig. (||) we see a 
reduction of 1/Z with increasing x. 
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From these results, obtained phenomenologically, we are able to conclude that 
Lai-^Sr^TiOa for x > 0.06 is a correlated metal which can be modeled with a local 
self-energy. The origin of the narrow quasiparticle peak and the lower Hubbard 
band can only be explained within a microscopic approach, e.g., DMFT. ThiiLis 
indeed possible, as has been shown recently using the LDA+DMFT approach.lia 

The Z-factor determined above is of the same order for both forms of the phe- 
nomenological self-energy [([?]) and fll2])]. Nevertheless, we find that the particular 
value of Z depends on the precise form of The dependence of Z on x also 

turns out to be different for the two self-energies. Using Eq. (^) its behavior is 
shown in the inset to Fig. ([!]). From Eq. ( |l2] ) we obtained values between 3.5 to 4 
for the corresponding x. This is because the phenomenological self-energy is used to 
fit the spectrum in a finite frequency window around ui = ^, which might be much 
larger then the actual frequency window for which Fermi liquid theory is valid. One 
should therefore be cautious in the interpretation of the actual values obtained for 
the effective mass m* . 

We also note an intriguing feature of the experimental data for x > 0.06.i 
When the curves are normalized to unity and plotted in a single figure, they are 
almost identical. In particular, the distance between the coherent peak and the 
top of the Hubbard band does not depend on x. This is in striking contrast to 
numerical calculations for a doped Hubbard model where doping inevitably leads 
to a shift of the chemical potential towards the lower Hubbard band (at least under 
the assumption of a homogeneous phase). Such a discrepancy has recently also 
been observed for the related compound Gdi-^Sr^TiOs where it was attributed 
to an inhomogeneous sample composition due a chemical phase separation into 
strongly and poorly doped domains £j Such a phase separation may be restricted 
to the surface, as bulk-sensitive experiments did not observe it. A possible difference 
between surface and bulk electronic structure in oxidic perovskites has indeed been 
reported for doped vanadates 

5. Conclusions and Final Remarks 

In this paper we presented a simple analytical expression (^|) for the self-energy 
which can be used to fit experimental data for strongly correlated electron systems. 
As an application we presented fits to the photoemission data for Lai-^Sr^TiOa, 
for which the phenomenological fit parameters all take reasonable values. 

There are certain difficulties in carrying out such a program, mainly due to the 
limited spectroscopic information which is currently available for correlated elec- 
tron systems. Direct and inverse photoemission data would help to enhance the 
quality of the fits, in particular, if they allow one to fit different bands. Angu- 
lar resolved photoemission experiments would give further information, and could 
provide insights into the validity of assuming a purely local self-energy. 

At the moment, our approach is particularly useful to fit data for systems with 
partially filled, degenerate bands which arc well separated from other completely 
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filled or empty bands as, e.g., in Lai-^Sr^TiOs. 
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